Enhancement of spin coherence in a spin-1 Bose condensate by dynamical decoupling 

approaches 



Bo- Yuan Ning, 1 Jun Zhuang, 1 J. Q. You, 2 and Wenxian Zhang^B 

1 Key Laboratory of Micro and Nana Photonic Structures (Ministry of Education), 
Department of Optical Science and Engineering, Fudan University, Shanghai 200433, China 
2 Department of Physics, Fudan University, Shanghai 200433, China 
(Dated: January 15, 2013) 

We study the enhancement of spin coherence with periodic, concatenated, or Uhrig dynamical 
decoupling iV-pulse sequences in a spin-1 Bose condensate, where the intrinsic dynamical instability 
in such a ferromagnetically interacting condensate causes spin decoherence and eventually leads to 
a multiple spatial-domain structure or a spin texture. Our results show that all the three sequences 
successfully enhance the spin coherence by pushing the wave vector of the most unstable mode in 
the condensate to a larger value. Among the three sequences with the same number of pulses, the 
concatenated one shows the best performance in preserving the spin coherence. More interestingly, 
we find that all the three sequences exactly follow the same enhancement law, fc-T 1 ^ 2 = c, with 
k- the wave vector of the most unstable mode, T the sequence period, and c a sequence-dependent 
constant. Such a law between fe_ and T is also derived analytically for an attractive scalar Bose 
condensate subjecting to a periodic dynamical decoupling sequence. 

PACS numbers: 03.75.Gg, 03.75.Kk, 03.67.Pp 



I. INTRODUCTION 

Spin coherence is of key importance to quantum de- 
vices utilizing spin degrees of freedom, e.g., a weak field 
magnetometer based on spin-1 Bosc-Einstein condensates 
(BECs), which requires the spin coherence time to be as 
long as possible in order to improve the device sensi- 
tivity [l|. However, due to the ferromagnetic spin ex- 
change interaction in the 87 Rb Bose condensates, the in- 
trinsic dynamical instability causes the unstable collec- 
tive modes to grow exponentially into a spin-domain or 
spin-texture structure and eventually sabotages the spin 
coherence of the condensates 0-0] • It is highly demanded 
to find a practical way to enhance the spin coherence 
in spinor Bose condensates in order to develop a higher 
precision magnetometer. In fact, the spin coherence is 
enhanced in a 8 7Rb spin-1 condensate by localizing the 
spin dynamics, which was proposed to unmask other in- 
teresting effects induced by the weak magnetic dipolar 
interactions pill]. 

Some experimental techniques may be employed to ma- 
nipulate the atomic spin dynamics of a spinor conden- 
sate. One is the magnetic radio-frequency pulses, which 
can rotate the spin in an arbitrary angle along a certain 
direction [10. Unfortunately, the spin exchange inter- 
action in spin-1 condensates is rotationally invariant so 
that the magnetic pulses are not able to control the spin- 
exchange-interaction induced effects. The other is the 
optical Feshbach resonance technique, which can adjust 
both the sign and the magnitude of the spin exchange 
coefficient ci = 4irti 2 (a2 — ao)/3m in a certain parame- 
ter range [lfj, where <zo,2 is the s-wave scattering length 
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and m the mass of the atom. Theoretically, it has been 
shown that the dynamical instability of the 87 Rb spin-1 
condensate can be suppressed by periodically changing 
the spin exchange coefficient in a cosine form [1]. 

From the viewpoint of quantum control theory, other 
forms of modulating the spin exchange coefficient c%(t) 
are not only easy to realize in experiments but also pow- 
erful to suppress the spin decoherence. Recently, Uhrig 
proposed an optimal sequence with N 7r-pulses to dy- 
namically decouple a qubit from its boson environment 
in a spin-boson model 11| . The Uhrig dynamical decou- 
pling (UDD) sequence has also been applied in a system 
of a center spin coupled with many nuclear spins [H, [HI] ■ 
The results [13j show that the UDD sequence can be used 
to remove the dephasing terms up to 0(T N+1 ), where 
T is the sequence period. This has aroused consider- 
able experimental interests to demonstrate the advan- 
tages of UDD (see, e.g., dMl). 

Moreover, the concate- 
nated dynamical decoupling (CDD) sequence has also 
been demonstrated to be powerful in suppressing the de- 
coherence in quantum systems [l7Hl9l] . Despite of the 
success of these dynamical decoupling (DD) sequences, 
they have never been tested in a Bose-Einstein conden- 
sate, where the decoherence is dominated by a very dif- 
ferent bath consisting of unstable collective modes. 

In this paper, we theoretically investigate periodic dy- 
namical decoupling (PDD), CDD, and UDD sequences 
on preserving the spin coherence of a homogenous 87 Rb 
spin-1 Bose condensate. By comparing the performance 
of these sequences, we identify the optimal sequence in 
suppressing the dynamical instability of the spin-1 Bose 
condensate. Our results show that all sequences sup- 
press the dynamic instability and the wave vector of 
the most unstable collective modes fc_ is pushed to a 
larger value. Although the CDD sequences behave bet- 
ter than the PDD or UDD sequences at the same num- 
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ber of pulses, the same law k-T 1 / 2 = c, with c being a 
sequence-dependent constant, is followed for all the three 
sequences. Moreover, in a scalar Bose condensate, our 
analytical result for a PDD sequence further verifies this 
law between fc_ and T. 

The paper is organized as follows. In Sec. [Ill we formu- 
late the spin dynamics of a spin-1 condensate subjecting 
to a dynamical decoupling pulse sequence and describe 
our metric to measure the performance of the sequence. 
In Sec. IIIIl we present numerical results for a spin-1 con- 
densate and compare the performances of PDD, UDD, 
and CDD sequences. To understand our numerical re- 
sults, we analytically investigate, in Sec. IIV1 a simple 
case of a scalar condensate subjecting to a PDD sequence. 
The conclusion is given in Sec. [V] and a brief derivation 
of our analytical results is provided in Appendix [A"l 



II. SPIN DYNAMICS OF A SPIN-1 
CONDENSATE SUBJECTING TO DYNAMICAL 
DECOUPLING SEQUENCES 

In the mean-field theory, the evolution of a homo- 
geneous three-component spin-1 Bose condensate is de- 
scribed by three coupled Gross-Pitaevski equations [20I — 
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where Hq = — (fi, 2 /2m)V 2 + V ex t + con is the spin- 
independent Hamiltonian. In the above equations, cq — 
4irh 2 (ao + 2a 2 )/(3m) represents the spin-independent in- 
teraction and nj(r) = |$i(r)| 2 is the total density of the 
condensate with i = 0,±1 and the condensate wave 
function of the ith component. 

The spin dynamics of the spin-1 condensate is deter- 
mined by the rest spin exchange terms containing c 2 . De- 
pending on the sign of c 2 , the spin exchange interaction 
can be either ferromagnetic (c 2 < 0, e.g., 87 Rb conden- 
sates) or antiferromagnetic (c2 > 0, e.g., 23 Na conden- 
sates) (23l - [3Cl | . It was shown that a ferromagnetically in- 
teracting spin-1 Bose condensate is dynamically unstable 
while an antiferromagnetically interacting one is dynam- 
ically stable Q. In other words, the spin coherence of 
a ferromagnetically interacting spin-1 condensate is lost, 
since some collective modes grow exponentially due to 
the intrinsic dynamical instability. 

According to the Bogoliubov approximation, the time 
dependent equations of motion for the collective modes in 
a condensate can be cast in a matrix form as follows [3l| : 
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ih- = M ■ x, 
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where x = (S^f +1 ,8^ , <S>E+i, 5^, S^*^) 11 denotes 

the amplitude of the collective modes and M is an effec- 



tive Hamiltonian. In general, M has a matrix form as 

M=[ I- (3) 



where for a spin-1 condensate 

A = s k I + c A + c 2 A 2 , 

B = c B + c 2 B 2 . 

Here, Sk = H 2 k 2 /(2m) is the kinetic energy of the collec- 
tive mode in a plane wave form with a wave vector k, I 
is a 3 x 3 identity matrix, and ^4o,2,-Bo,2 are also 3x3 
matrices given by 
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Obviously, matrix M is neither Hermitian nor symmet- 
ric although Ao, 2 are Hermitian and i?o,2 are symmetric. 
However, the matrix M has paired eigenvalues with the 
same amplitude but different sign, due to its particular 
form 31|. The dynamical instability of a condensate is 
manifested as the existence of at least one pair of complex 
eigenvalues [|[. 

By employing the optical Feshbach resonance tech- 
nique, the spin exchange coefficient c 2 can be tuned not 
only in amplitude but also in sign It occurs natu- 
rally to utilize such a technique to suppress the dynamical 
instability in a ferromagnetically interacting spin-1 Bose 
condensate. Previous theoretical investigation has indeed 
demonstrated the suppression effect by modulating c 2 (t) 
in a simple cosine form [8|. 

From the viewpoint of the control theory, there are 
many more complicated DD pulse sequences which may 
outperform the simple cosine modulation of c 2 (t) in sup- 
pressing the dynamical instability and enhancing the spin 
coherence. In addition, these pulse sequences are easier 
to implement in experiments. Hereafter, applying a dy- 
namical decoupling pulse means alternating the sign of 
c 2 but keeping its amplitude unchanged. 

A PDD sequence denotes flipping the sign of c 2 back 
and forth periodically [see top panel in Fig. [Tj . By con- 
catenating PDD sequences at different levels, we obtain 
the CDD sequences [see middle panel in Fig.[Tj. A UDD 
sequence [see bottom panel in Fig. [T] containing N pulses 
is defined as EMI 
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FIG. 1. Various PDD (top), CDD (middle), and UDD (bot- 
tom) sequences. The vertical short lines mark the time of 
each applied pulse. 

where Tj is the delay between the (j — l)th and jth pulses 
and T — ^2^1 Tj is the period of the total sequence. 
We notice that the CDDi (CDD at one-level concatena- 
tion), UDDi (one-pulse UDD), and PDD sequence are 
the same, and CDD2 and UDD2 are also the same. 

There are various choices of metric to measure the 
pulse sequence effect. We choose the wave vector of the 
most unstable mode, fc_, since this mode dominates in 
the dynamics of the spin-1 condensate The larger the 
value of fc_ is, the better the performance of the pulse 
sequence will be. For a spin-1 condensate under many 
DD pulses, we adopt a numerical method (se below) to 
obtain the fc_ . 

By applying a sequence of N DD pulses, the time evo- 
lution operator Ut of the collective modes of the spin-1 
Bose condensate is 



Ut = Texp 



dt'M{t') 



— exp[— iM±TN + i] . . . exp[— iM^T%\ exp[— iM±T\], 

(5) 

where T denotes the time-ordering operator, and M± 
denotes the M matrix with C2 being either positive or 
negative. By diagonalizing Ut and reexpressing it as 
Ut = exp(—iTF)V with F a 6 x 6 diagonal matrix, 
we identify the wave vector of the most unstable mode fc_ 
by searching the largest imaginary part of F for different 
wave vector k M. 



III. NUMERICAL RESULTS FOR A SPIN-1 
BOSE CONDENSATE 

Following the procedure described at the end of Sec. ITU 
we numerically obtain the dependence of the unstable 
modes on the DD sequence periods for a spin-1 Bose con- 
densate. For the numerical simplicity, we set m = 1 and 
H = 1. The results are presented in Fig. [5] and Fig. 
Several interesting observations are listed as follows: 
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FIG. 2. (Color online) Dependence of the imaginary part of 
the Floquet operator F on the wave vector k and the pe- 
riod T for sequences (a) PDD 3 , (b) UDD 5 , and (c) CDD 3 . 
Lower panels show the k dependence for (d) PDD3, (e) 
UDD 5 , and (f) CDD 3 at T = T /20 which is marked by 
the vertical white lines in panel (a), (b), and (c), respec- 
tively. Here, we use parameters [8j, [34] ko = 1 and To = 

7r/ C2n-y/no(i — no) cos(0o/2) = 7r/(0.4|c2 \n), where the ini- 
tial condition is 60 — and no = 0.8. In calculations we have 
set Co = 1, n = 1, and C2 = —0.0046 which are the parameters 
of 87 Rb spin-1 condensates. 



i) Compared to the free evolution (T = To), all the 
DD sequences suppress the dynamical instability of 
the spin-1 condensate in three ways, narrowing the 
area of the unstable region, shifting fc_ to a larger 
value, and decreasing the value of F_ = Im(F)\ k_ 
(the exponential growth rate of the most unstable 
mode). 

ii) Comparing the performances of the same DD se- 
quence at different period T, i 7 ! is almost the same 
but k- is shifted to a larger value when decreasing 
T. 

iii) At the same number of pulses (N — 5) and the 
same sequence period, the UDD sequence shifts fc_ 
the least, while the CDD and PDD sequences be- 
have almost the same (but both better than the 
UDD) [see Fig.^a)]. Moreover the CDD sequence 
suppresses the growth rate F_ more than the PDD 
[see Figs.Hd) and df)]. 

iv) The relationship between fc_ and T for all PDD, 
UDD, and CDD sequences follows exactly the same 
power law: k-T 1 / 2 = c, where c is a sequence- 
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dependent constant and can be obtained by numer- 
ically fitting the slope of the lines in Fig. 03a). 

The above results clearly indicate that the dynamical 
instability in a ferromagnetically interacting spin-1 Bosc 
condensate is essentially inhibited once the wavelength 
of the most unstable mode 27r/fc_ is less than the spin 
healing length £ = hj ^/2m\ci\n. Our results shows that 
CDD sequences most effectively suppress the dynamical 
instability and they are the best choice to enhance the 
spin coherence in a spin-1 Bose condensate. 

It is not easy to fully understand all the above nu- 
merical results for a spin-1 Bose condensate in a simple 
picture, especially the result in (iv), because it is diffi- 
cult to obtain analytical results with either the effective 
Hamiltonian M matrix [Eq. or the F matrix. To 
circumvent these difficulties, we investigate the effect of 
the same DD sequences but in a scalar Bose condensate, 
where it becomes possible for us to obtain the same be- 
haviors but in an analytical way. 



results for a PDDi sequence. For convenience, we reform 
the matrix M± in Pauli matrices as 



M± = (s k ± gn) ■ a z ± ign ■ a v . 



(7) 



By utilizing the properties of Pauli matrices, it is 
straightforward to obtain the time evolution operator 

U± = exp(— itM±) 

= cos(tx±) [(eic ± gn) ■ a z ± ign ■ a y ] sin(rx±), 

x± 



where x± — y (e^ ± 2gn) and r = T/2. The propaga- 
tor in a complete period is 

U * = U - U + = (-^+iQ v ~h + iQ, V )> (8) 



and the corresponding eigenvalues of Ut is 



P±JQl + Ql-Ql 



(9) 



where 



IV. RESULTS FOR SCALAR BOSE 
CONDENSATE 

Similar to the case of a spin-1 Bose condensate, we 
also focus on the change of fc_ for the most unstable 
mode of a scalar condensate subjecting to DD pulses, 
which alternates the sign of the interaction coefficient 
g = Airti?a s /m, with a s the s-wave scattering length. 
The interaction between atoms is repulsive if g > and 
attractive if g < [33]. For a homogeneous scalar con- 
densate with attractive interaction, the system is intrinsi- 
cally unstable [34| . By periodically flipping the sign of g, 
it is possible to suppress the instability of an attractively 
interacting Bose condensate [gj]. Fortunately, both the 
sign and the amplitude of g can be precisely tuned in ex- 
periments via the Feshbach resonance technique (35l . |36| . 

The time-dependent effective Hamiltonian M (t) ma- 
trix of the collective modes in the scalar condensate mod- 
ulated by a PDDi sequence is given by 



M(t) 



where 
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M 4 
M_ 



£fc ± gn 
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t e 
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Following exactly the same procedures in Eq. ([5]) for the 
above M matrix, we first diagonalize the corresponding 
time evolution operator Ut, then find out the eigenval- 
ues of the Flouqet operator F, and finally extract the 
mode containing largest imaginary part F— and the cor- 
responding wave vector k_ . For a scalar Bose condensate, 
the matrix M(t) is much simpler than the one in a spin-1 
condensate, so it is possible for us to obtain analytical 



P = cos(rx_) cos(ra;+) — sin(ra;_) sin(ra;+). 

X-X+ 
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X-X+ 

^ 9 n ■ , \ i \ 9 n i \ ■ i \ 
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r. ^k- gn 
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Since A = exp(— iTF), the imaginary part of the Flou- 
quet operator follows as 



1 



Jm(F) = - In |A|. 



(10) 



We notice that the kinetic energy e k of the most un- 
stable mode is much larger than the interaction energy 
gn if the PDDx sequence period T is small. By assum- 
ing gn -c Sk and keeping to the second order of gn/sk, 
we analytically obtain the relationship between the wave 
vector fc_ of the most unstable mode and the sequence 
period T from Eqs. ©- (HSJ as 



c, 



(11) 



with c = y/2n. For a detailed derivation of Eq. ([TT|) , see 
Appendix [X] 

It is difficult to obtain the analytical results for other 
more complex sequences such as UDD5, PDD3 and 
CDD3, so we carry out numerical calculations for the 
scalar condensates. The results are shown in Fig. EJb). 
It is clearly shown in Fig.[3]that no matter in a spin-1 con- 
densate or in a scalar condensate, all sequences exhibit 
the same relationship between k- and T: k^T 1 ! 2 = c, 
with c a sequence-dependent constant (see Table A} . 
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FIG. 3. (Color online) Numerical results of power law rela- 
tionship between the wave vector k- of the most unstable 
mode and the sequence period T for (a) spin-1 condensates 
and (b) scalar condensates: open squares — CDD3, crosses — 
PDD3, filled circles — UDD 5 , filled triangles — PDDi. The 
specific fitting results are shown in Table [U The parameters 
for the scalar condensates are g — —0.0046 and n = 1. 



TABLE I. Numerical fitting values of the slope, intercept and 
calculated c 2 /2tv for lines in Fig.[3]for spin-1 condensates [de- 
noted below as (a)] and scalar condensates [denoted as (b)]. In 

calculations, c 2 = 7rexp(2yo)/ \\c2 \n\J no(l — no)J for spin-1 

condensates and c 2 = 2tt &cp(2yo) / {\g\n) for scalar conden- 
sates. 





slope 


intercept (yo) 


c 2 /2tt 


(a) PDDi 


-0.473 


-3.051 


0.608 


(a) UDD 5 


-0.499 


-2.799 


1.006 


(a) PDD 3 


-0.495 


-2.583 


1.551 


(a) CDD3 


-0.495 


-2.583 


1.551 


(b) PDDi 


-0.498 


-2.681 


1.020 


(b) UDD 5 


-0.499 


-2.341 


2.013 


(b) PDD 3 


-0.499 


-2.137 


3.027 


(b) CDD 3 


-0.499 


-2.137 


3.027 



the unstable collective modes in a spin-1 Bose Einstein 
condensate. Our numerical results show that the three 
DD sequences can be used to effectively protect the co- 
herence of the spin-1 condensate. The wave vector fc_ 
of the most unstable mode is shifted to a much larger 
value, which in fact enhances the spin coherence of the 
condensate. Furthermore, using the analytical results in 
a scalar condensate, we reveal an interesting relation be- 
tween k- and sequence period T: k^T 1 / 2 — c, where c is 
a sequence-dependent constant. 
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Appendix A: Derivation of Eq. (|T1[) 



According to the discussion in Sec. IIV1 we extract the 
most unstable mode which has the largest growth rates 
F_, satisfying d[Im(F)]/dk\ k _ = 0. Using Eq. (JlO]), we 
find that the derivative of F is 



dlm(F) _ 1 d(ln|A| 



dk 



T dk 



(Al) 



where 



|A| = 



P 2 + Ql-Ql-Q 



A is real, 

A is complex, 



(A2) 



V. CONCLUSION 

In conclusion, we have theoretically investigated the 
PDD, CDD and UDD pulse sequences on suppressing 



with P and Q x ,y,z given in Eq. @. In the complex A 
case, |A| = 1 and Im(F) — 0, which means a stable 
collective mode. While in the real A case, by defining 
e = gn/sk 1, A is reformed as 



6 



A = cos 2/+ cos y_ ^=^= sin y + sin w_ 

VI - 4e 2 



k 2 -2 



1 — cos 2y + cos 2y_ + 2e 2 (cos 2y + + cos 2y_ + cos 2y + cos 2y_ — 3) + yl — 4e 2 sin 2y + sin 2j/_ 

'(A3) 

where ?/± = efcT-^/l ± 2e. By keeping terms up to e 2 , we obtain 



A R3 cos(y + + y_) — 2e 2 siny + siny_ + *J— sin 2 (y + + y_) — 4e 2 cos(y + + y_) siny + sin?/_ . (A4) 



Considering the fact that e 1 and A is real (which Another concern is the suppression of F- . According 

requires the terms under the square root to be of the to Eq. (|A4[) . when k-T 1 / 2 = ^/2tt, 
order of e 2 ), we have y + ~ y_ ~ k 2 r/2, and 

x < jl \.nf\ fAK\ A = -l + 4^ (A6) 

Ass cos(y+ +0(e). (A5) fc 2 v ; 

Thus, as long as y + + y_ = 7r, if we keep terms up to the order of e. Finally, we obtain 

that 

dFL 

dk ' P -I^h^M* 



F_ = -ln^l-4^, (A7) 
which gives rise to the relationship in Eq. i.e., 

k_T x ^ 2 — \/2tt. which shows that .FL depends very weakly on k-. 
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